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BIRTH OF LIMIT CYCLES FOR A CLASS OF 
CONTINUOUS AND DISCONTINUOUS 
DIFFERENTIAL SYSTEMS IN (d + 2)-DIMENSION 

JAUME LLIBRE * 1 , MARCO A. TEIXEIRA 2 AND IRIS O. ZELI 2 


Abstract. The orbits of the reversible differential system x = 
—y, y = x, z = 0, with x,y € R and z € R d , are periodic with 
the exception of the equilibrium points (0,0, z). We compute the 
maximum number of limit cycles which bifurcate from the periodic 
orbits of the system x = —y, y = x, z = 0, using the averaging 
theory of first order, when this system is perturbed, first inside 
the class of all polynomial differential systems of degree n, and 
second inside the class of all discontinuous piecewise polynomial 
differential systems of degree n with two pieces, one in y > 0 and 
the other in y < 0. In the first case this maximum number is 
n d (n — l)/2, and in the second is n d+1 . 


1. Introduction and statements of the main results 

Limit cycles have been used to model the behavior of many real pro¬ 
cess and different modern devices. In general to prove the existence 
of limit cycles is a very difficult problem. One way to produce limit 
cycles is perturbing differential systems that have a linear center. In 
this case, the limit cycles in a perturbed system bifurcate from the pe¬ 
riodic orbits of the unperturbed center. The search for the maximum 
number of limit cycles that polynomial differential systems of a given 
degree can have is part of 10 th Hilbert’s Problem and many contribu¬ 
tions have been made in this direction, see for instance HU US EE] and 
the references quoted therein. 

Recently the theory of limit cycles has also been studied in discon¬ 
tinuous piecewise differential systems. The analysis of these systems 
can be traced from Andronov et al. [1] and still continues to receive 
attention by researchers. Discontinuous piecewise differential systems 
is a subject that have been developed very fast due to its strong appli¬ 
cations to other branches of science. Currently such systems are one of 
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the connections between mathematics, physics and engineering. These 
systems model several phenomena in control systems, impact in me¬ 
chanical systems, nonlinear oscillations and economics see for instance 
EEUEBj • Recently they have been shown to be also relevant as 
idealized models for biology [15] and models of cell activity [7 : , f32l 133] . 
For more details see Teixeira [3TJ and all references therein. 

As we have said it is not simple to determine the existence of limit 
cycles in a differential system. The simplest case for determining limit 
cycles is in planar continuous piecewise linear systems when they have 
only two linear differential systems separated by a straight line. Even 
in this simple case, only after a delicate analysis it was possible to show 
the existence of at most of one limit cycle for such systems, see [8j or 
an easier proof in [20]. 

Planar discontinuous piecewise linear differential systems with only 
two linearity regions separated by a straight line have been studied re¬ 
cently in [ITIT1 [T23 5 among other papers. In ra some results about the 
existence of two limit cycles appeared, so that the authors conjectured 
that the maximum number of limit cycles for this class of piecewise 
linear differential systems is exactly two. However in [12] strong nu¬ 
merical evidence about the existence of three limit cycles was obtained. 
As far we know the example in [T2] represents the first discontinuous 
piecewise linear differential system with two zones with 3 limit cycles 
surrounding a unique equilibrium. Recently in [22] it is proved that 
such a system really has three limit cycles. 

There are several papers studying the limit cycles of the continuous 
piecewise linear differential systems in M 3 , see for instance [3I2DI22IMI 
[25] . Our goal is study the periodic solutions of discontinuous piecewise 
polynomial differential systems in More precisely the objective 

of this paper is to study the existence of limit cycles in continuous 
and discontinuous piecewise polynomial differential systems in M d+2 , 
where the discontinuous differential system has two zones of continuity 
separated by a hyperplane. Without loss of generality we shall assume 
that the set of discontinuity is the hyperplane y — 0 in M d+2 . So we 
consider the linear differential system in R d+2 given by 


x = -y 

(1) y = x 

ze = 0 

for £ = 1,... ,d and x,y G M, z e M d , where the dot denotes de¬ 
rivative with respect to the time t , which is reversible with respect to 
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4>(x,y,z) = ( x,—y,z ) . First we are interested in studying the exis¬ 
tence of limit cycles of the continuous polynomial differential system 
given by 

x = - y + eP a (x,y,z), 

( 2 ) y = x + eP b {x,y,z), 

zi = eP Ce (x,y,z), 

and after we also study the existence of limit cycles of the discontinu¬ 
ous piecewise polynomial differential system formed by two polynomial 
differential systems separated by the hyperplane y = 0 , namely 


(3) 


X = 

-y + eP a (x,y,z), 

y = 

x + eP b (x,y,z), 

Z£ = 

eP Ce {x,y,z), 

X = 

-y + eQ a (x,y,z), 

y = 

x + eQ b (x,y,z), 

z e = 

£Qct(x, y, z), 


if y > 0 


if y < 0 


where £ 7 ^ 0 is a small parameter and t=\,,d. I 11 this systems the 
polynomials P a , P b , P cv Q a , Qp, Q le are of degree n in the variables 
x, y and z, more precisely 


n n 


Pa{x,y,z ) 

= E 

a ijk x l y J z k , 

p b 1 

[x,y,z) = 

= ^2 b ijk x l y J z k , 


i+j+k= 

=0 



i+j+k =0 

P ct (x,y,z ) 

= E 

cujkX 1 y° z k , 

Qc 

i(x,y,z) 

= ^2 a ijkX l y J Z k , 


i+j+k= 

=0 



i+j+k =0 

Qp(x,y,z) 

TL 

= E 

Pijk'X 1 y 3 Z k , 

Q^ e (x,y,z) 

Tl 

= '22 r y^jkX l y J z k 


i+j+k =0 i+j-\-k =0 


In this expressions k is a multi-index and i +j + k denotes i +j + k± + 
... + kd, z k denotes the product z^ 1 ... z kd with z — (zi,... ,za) G M d , 
and ciijk denotes the coefficient a'ijk 1 ...k d of x l y 3 z\ l ... z^ d . 

It is clear that systems ([2]) and (J3]) coincide for £ = 0 and they have 
linear centers at every plane z = constant. In this paper we establish 
for e 7 ^ 0 sufficiently small the maximum number of limit cycle of these 
systems that bifurcate from the periodic orbits of these linear centers 
using the averaging theory of first order. The following result presents 
the results for the continuous case. 
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Theorem 1. Using the averaging theory of first order for |e| 7 ^ 0 suffi¬ 
ciently small the maximum number of limit cycles of the polynomial dif¬ 
ferential system (J 2 ]) is at most n d {n — 1 )/ 2 , and this number is reached. 

In the next theorem we present results for the discontinuous piecewise 
polynomial differential system ( 1 H 1 ) . 

Theorem 2. Using the averaging theory of first order for |e| 7 ^ 0 suffi¬ 
ciently small the maximum number of limit cycles of the discontinuous 
piecewise polynomial differential system ((2D is at most n d+1 , and this 
number is reached. 

Corollary 3. Under the assumptions of Theorem Q] if additionally 
a 00 k = b 00k = 0 for all k, the limit cycles can be chosen as close to 
the origin of R d+2 as we want. 

Corollary 4. Under the assumptions of Theorem (K if additionally 
aook = b 00k = aook = fiook = 0 for all k, the number of limit cycles 
of the discontinuous piecewise polynomial differential system ( 131 ) is at 
most n d (n — 1), and this number is reached. Furthermore, the limit 
cycles can be chosen as close to the origin ofM. d+2 as we want. 

Corollaries [3] and Q] provides information of the Hopf bifurcation of 
systems (j2J) and Q . More precisely, Corollaries [3] and Q] show that at 
least n d (n — l )/2 and n d {n — 1) limit cycles of systems ((2D and ((2D can 
bifurcate from the origin of R d+2 , respectively. The results of Corollary 
[3] in the particular case n = 2 coincides with the result obtained in 
Theorem 1 of [26] . 

To prove these results we use the classical averaging theory, see for 
instance [30, (21] for a general introduction to this subject. This theory 
have been used for years to deal with continuous differential systems. 
The principle of averaging has been extended in many directions and 
recently in [TOJ the authors extend the averaging theory for detecting 
limit cycles of certain discontinuous piecewise differential systems, via 
the Brouwer degree and the regularization theory. 

As far as we know this method is one of the best methods for deter¬ 
mining limit cycles in discontinuous piecewise differential systems and 
has already been used by some authors. I 11 id the method is used 
for determining the maximum number of limit cycles that bifurcate 
from the periodic solutions of some family of isochronous cubic poly¬ 
nomial centers perturbed by discontinuous piecewise cubic polynomial 
differential systems with two zones separated by a straight line. In [18] 
limit cycles for discontinuous piecewise quadratic differential systems 
with two zones was studied using the averaging theory. Also in [29] 
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the averaging theory was applied to provide sufficient conditions for 
the existence of limit cycles of discontinuous perturbed planar centers 
when the discontinuity set is a union of regular curves. 

We have organized this paper as follows. In Section [2] we briefly 
present notation and basic concepts of the averaging theory of first or¬ 
der for continuous differential systems (see Theorem [5]) and for discon¬ 
tinuous differential systems (see Theorem 0 . In Section [3] we present 
the proof of results of this paper. More precisely, the proof of Theorem 
|T| is presented in the Subsection 13. 11 and in the Subsection 13. 21 we prove 
the Theorem [2l Finally in the Subsection 13.31 we presented the proofs 
of the Corollaries [3] and HJ 


2. Basic results in averaging theory 


In this section we present the basic results from the averaging theory 
of first order that we shall use for proving the results of this paper. 
The following theorem provides a method for studying the existence 
of periodic orbits of a differential system. For more details on the 
averaging method see for instance | 33 ]. 

Let D be an open subset of M n . We denote the points of M x D as 
and we take the variable t as the time. 

Theorem 5. Consider the differential system 


x = eF(t, x) + e 2 R[t, x, s), x(0) = 0, 


( 4 ) 


where F : K. x D —> M n and R : R x U x (— e, e) —> M n are continuous 
functions and T—periodic in the first variable. Define the averaging 
function f : D -A M n as 



and assume that 

(i) the functions F, R , D X F , D 2 X F and D X R are defined, con¬ 
tinuous and bounded by a constant M (independent of e) in 
[0, oo) x D and for e G (0, £o]> 

(ii) for p G D with f(p) = 0 we have \Jf(p)\ f 0, where \Jf(p)\ 
denotes the determinant of the Jacobian matrix of f evaluated 
at p. 

Then for |e| > 0 sufficiently small there exists a T—periodic solution 
x(t, e) of the system (J3j) such that x(0, e) —» p as e ^ 0. 

Now let h : Ixd —» M be a C 1 function with 0 6 I as a regular value, 
and E = h -1 (0). Let X, Y : M x D — » M n be two continuous vector 
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fields and assume that h, X and Y are T— periodic in the variable t. 
We define a discontinuous piecewise differential system as 


( 6 ) 


f X(t,x) if h(t,x)> 0, 
[ Y(t, x) if h(t, x) < 0. 


We rewrite the discontinuous differential system as follows. Consider 
the sign function defined in M\ { 0 } as 


sign(w) 


1 if u > 0 , 
-1 if u < 0 . 


Then system (J 6 ]) can be written as 


x — Z(t, x) — Fi(t,x) + sign ( h(t,x )) F 2 (t,x), 


where 

Fi{t,x) = ^ (X(t,x) + Y(t,x)) and F 2 (t,x) = ^ (X(t,x) - Y(t,x )). 

The following theorem is a version of Theorem [5] for studying the 
periodic solutions of discontinuous differential systems. 


Theorem 6 . Consider the discontinuous differential system 
(7) x = eF(t, x) + e 2 R(t, x, e), 

with 


F(t,x) = Fi(t, x) + sign (, h(t,x )) F 2 {t,x), 

R(t, x, e) = Ri(t, x, e) + sign (h(t, x)) R 2 (t, x, e), 

where Fi,F 2 : M x D — » M n , R\,R 2 : M x D x (— e,e) W 1 and 
/t:Ixh->R are continuous functions, T—periodic in the variable t. 
We also suppose that h is a C 1 function with 0 as a regular value and 
we denote E = h -1 (0). Define the averaged function f : D —* R" as 


( 8 ) 

and assume that 


f(x) = [ F(s,x ) ds, 
Jo 


(i) the functions Fi,F 2 ,Ri, R 2 and h are locally Lipschitz with re¬ 
spect to x; 

dh 

(ii) -Qjr{t,x) 0 for all (t,x) E E; 

(iii) for p E C with f(p) = 0, there exist a neighborhood U C C of p 
such that f(z) 0 for all z E U\ {p} and d B (/, U, 0) 7 ^ 0 (ds 
is the Brouwer degree of f in p). 
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Then for |e| > 0 sufficiently small there exists a T—periodic solutions 
x(t, x) of system (J7D such that x(t, e) —> p as e —* 0 . 

For a proof of Theorem [ 6 ] see Theorem A and Proposition 2 in 
[19]. Here we emphasize that if / in (J 8 J) is C 1 then the hypotheses 
ds{f, U, 0) f 0 holds if J f (p) \ 0, see for more details [27], 


3. Proof of the main results 

We devoted this section for the proof of results of this paper. For 
this we consider the notations introduced in the Introduction. 


3.1. Proof of Theorem [IJ Applying the change of variables 

(9) x = rcos9, y — r sin 9 and z = z, 
system (J2J) becomes 

n 

1 + - ^ r' l+j z k [aijk cos* 9 siid +1 9 + bij k cos * +1 9 sin - 7 9), 

i+j+k =0 
n 

e ^ r l+j z k ( atjk cos * +1 9 sird 9 + bij k cos* 9 sin J+1 9 ) , 

i+j+k =0 
n 

e ^ cujk r l+j z k cos* 9 sin - 7 9, 

i+j+k =0 

Essentially for e 0 sufficiently small we study the existence of limit 
cycles bifurcating from periodic orbits of this system when £ = 0, and 
that are contained in the cylindrical annulus 

(10) A — {(9,r, z) : r 0 < r < ri, 0 6 S 1 , z G W 1 }. 

So for £ sufficiently small 9 > 0 for every ( 9 , r, z) G A. 

Now taken as new independent variable 9 instead of t so f! 2 |l in A 
can be written as 


9 = 

r = 

k = 


r = 


( 11 ) 


£F 1 (9,r,z) + 0(£ 2 ) 
£F t+ i(9,r,z ) + 0(£ 2 ), 
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for £ = 1 ,d, where the prime denotes derivative with respect to 
the variable 9 , and 

Fi{0,r,z) 

( 12 ) 

F t+1 (d,r,z) 

for £ — 1, 2 ,..., d. 

System (fill) is 27r-periodic with respect to the independent variable 
9, and satisfies the hypotheses of the Theorem [5] for e 0 small and D 
fixed. Now we compute the averaged function / = (/i, / 2 , ■ ■ ■, fd+i) 
given in ([5]) with T = ‘Itx. 

Now taking 

H{p, q ) = / cos p 9 sin 9 9 d 9, 

Jo 


n 

= ^ r l+J z k (a i]k cos* +1 9 sird 9 

i+j+k =0 

+bijk cos* 9 sin 9+1 9 ), 

n 

= ^ r l+j z k cujk cos* 9 sin 9 9 , 

i+j+k =0 


note that H(p, q ) ^ 0 if only if p and q are simultaneously even. 


So we obtain 
(13) 

fi(r,z) 


fe+i(r,z) = 


r-2lt 


Fi(9, r, z) d 9 


n 


n 

r l+3 z k a ijk 

i + j + k = 0 
i odd , j even 


+ y r *+l z k b ijk fj, {iJ+ 1 ) 


i + j + k = 0 
i even , j odd 


-2it 


F e+ i(9,r, z) d0 


n 

^ ) T 3 Z Cgij k 

i + j + k = 0 
even 


for i — 1 , 2 ,..., d. 

We split the proof in two parts. First assume n is odd, so 


fi(r, z) = Apr + A 3 r 3 + ... + A n r n 


where 

(14) 
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n—p 


a p = J2 z " 

k =0 




^ ^ ®ijk /-ffi+l ,j) 4“ ^ ) bijk h(z.j+l) 


i + j = P 
\i odd,j even 


i+ j =P 
i even,j odd 


7 


for p = 1 , 2 ,..., n. So we write fi = rj\ with 

fi(r, z)= Ai + A 3 r 2 + ... + A n r n ~ l . 


Since r > 0 it is sufficient to solve (/i, / 2 ,..., fd+i) = (0,..., 0) to 
determine the number of solutions of / = 0 . As f\ is a polynomial in 
the variables r and z G of degree n — 1 and fi + ± are polynomials in 
the variables r and z 6 R d of degree n for l — 1, 2,..., d, by Bezout’s 
theorem (see i) (/i , / 2 ,..., fd+i) has at most n d (n — 1 ) solutions. 
However fi is even on variable r then we consider only solutions with 
r > 0 , in this case the maximum number of solutions of / = 0 is 
n d (n — l)/ 2 . 

Now we prove that this number is reached. For this, we exhibit a 
particular case for which this occurs. Let 0 , bij 0 ^ 0 and we take zero 
all the other ciijk,bijk, then /i(r, z) is a real polynomial that does not 
depend of z E of degree n — 1, that is 


fi^A 1 + A 3 r 2 + ... + A n r n ~ x 


where 


A p — 


a v° d(i+l,j) + 


E 


bijO d(i,j +!)• 


i+j=p 
i odd, j even 


i + j=p 
i even, j odd 


On the other hand, we take cajk = 0 if j, k\, k t _ i, ki + i, ..., k d ^ 0 
for each i — 1 , 2 ,..., d so that 

n 

fi+i(r,z) = y zt l (2vr c mkl ). 

ke =0 


In this particular case, we can take a^o, bij 0 and Q 0 o^ hi order that all 
coefficients of f\ and fi + 1 are linearly independent for all i = 1,2 ,,d. 
So we can choose these coefficients in such a way that /i has (n — l )/2 
simple positive real roots and fe + ± has n simple real roots for each 
i = 1, 2,..., d. Then (/], / 2 ,..., fd+i) has n d (n — l)/2 solutions with 
r > 0 . 
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Now assume n is even. Then / = rfi where 

fi(r, z) = A 1 + A 3 r 2 + ... + A n _ir n ~ 2 . 
with A p given in (fT4lh 

Note that fi has degree n — 1 as a polynomial in the variables r and 
z and ff + i given in (TT3]l has degree n as polynomials in the variables r 
and z for all l — 1, 2 , ..., d, so (/i, / 2 ,..., fd+i) have at most n d {n—1)/2 
solutions of type (r, z) with r > 0. 

To prove that this number is reached we consider a\Q kl , &oifci 7^ 0 and 
we take zero all the other a^, %fc- Then 

n —1 

fi( r ,z) — Ai — 'Y / z^ 1 (aiofci ^( 2 , 0 ) + feoifci AKo, 2 )) 

fci=0 

is a complete polynomial on variable zi of degree rt — 1. Now for 
£ = 2, 3..., d we take ci^o, Qoofc £ 7^ 0 and we take zero all the other 
C£ijk so that 

n 

f 2 (r,z) = rl+J Cl u° 

i + j = 0 

i,j even 
n 

ft+i(r,z ) = Y, zY c mke 2n, 

k e =0 

for l — 2, 3 ..., d. 

Here we take a wkl , b 0 i kl , cuj 0 and Cmk e in such a way for that 
all coefficients of f± and fe + ± for t = 1,2 ,,d, are linearly inde¬ 
pendent. Therefore we can choose these coefficients in order that fi 
has 7i—l simple real roots, / 2 has n/2 simple positive real roots and 
j f+ j has n simple real roots for each £ — 2, 3,..., d . In this case 
(/i, / 2 ,..., fd+i) = (0,..., 0) has n d (n — l)/2 solutions with r > 0. 
Furthermore by independence of the coefficients these solutions can be 
taken in a way that the Jacobian of / in all these solutions is nonzero. 
This completes the proof of Theorem [0 

3.2. Proof of Theorem [2j Analogously to what we did in the previ¬ 
ous section we shall study the existence of limit cycles bifurcating of 
periodic solutions of system ((2D when £ = 0 which are contained in the 
cylindrical annulus A defined in (TTCTT) . 

In A we have for e sufficiently small 9 > 0 for all (0, r, z) G A. Taking 
9 as independent variable system ((3j) in A becomes 
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( 15 ) 


r 

z 


£F 1 (9,r,z) + 0(£ 2 ), j 
£F £+1 (9,r,z) + 0(£ 2 ) j 

sGt(0, r, z) + 0(£ 2 ), 1 

£Gi+\{0i r, z) + 0(£ 2 ) J 


if h(9, r, 2 ) > 0, 
if h{9 , r, z) < 0, 


where 


(16) 


n 

Gi ( 9 , r, 2) = r i +j z k (ctijk cos * +1 9 sin J 6 * 

*+j+fc=0 

+/3ijk cos* 9 sin J+1 9 ), 

n 

G e+ 1 (#, r, 0) = ^ r l+ 3 z k ^ujk cos* 6 1 sin J #, 

i+j+k =0 


for £ = 1, 2,..., d, F] and F 2 given in (1T2|) . and h(9, r, z) 

dh 

Then for ( 9,r,z ) E h _1 (0) we have —(9,r, z)\g e { 0)7T y = cos# 


±1 ^ 0 . 


d9 


= sin#. 

|0€{O,7r} = 


Now, let 


/*7T 

I{ p ,q) = / cos p # sin 9 # d# and J( P , q ) 
Jo 


"2tt 


cos p 9 sin 9 # d#. 
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Then (—1 ) 9 J( p , q ) = I(p,q) and I( p , q ) = J( P , q ) = 0 if only if p is odd. 
Thus we have 

/ it p2ir 

F 1 (6,r,z)d6 + Gi(0, r, z) d6 


^ r l+J z k (a ijk + (-1 ) 3 a ijk ) J (i+ 


i + j + k = 0 
i odd 


(17) 


+ T r' + ’z k ( b ijt - (-1 YPijk) /(u + i), 


i + j + k = 0 
i even 


f* 2n 


fe+i(r,z)= F e +i(6,r,z) d6+ G e+1 (9,r,z) dd 


y r l+J z k (ctijk + (-1 yiujk) I(i,j ), 


i + j + k = 0 
i even 


for £ — 1, 2,..., d. 

Then fg + 1 are polynomials on variables r and z of degree n for each 
£ — 0,1,..,, d. By Bezout’s theorem / has at most n d+1 solutions (with 
r > 0). 

To prove that this number is reached we choose a particular example. 
So take a^o — (—l)-'ayo ) hjo — (—l) J Ajo 7^ 0 and we take zero all the 
other aij k ~ (—1 YcXijk, hjk ~ (—1 Then f\ is a real polynomial on 

r of degree n that does not depend of z. 

Analogously for each £ — 1, 2,..., d we take Q 0 ofc f + 7ioofc« 7^ 0 and 
we take zero all the other c^jk + 'yujk, so that 

n 

fe+i( r -> z ) — Zj e (cmk e + r Ymk e )'K- 

ke =0 

Under such conditions all coefficients of f\ and ft + 1 for £ = 1, 2,..., d 
can be taken to be linearly independent from the appropriate choice 
of a.ij o, a ij0 , b ij0 , Ajo, Qoo^ and ~/mok r Such values can be taken so 
that f\ is a complete real polynomial on variable r of degree n, and 
therefore it can have n simple positive real roots. Furthermore for each 
l = 1,2 ,,d the polynomials fg + 1 is a complete real polynomial on 
variable zg with n simple real roots. So for this particular case the 
number of zeros of / is n d+1 . By independence of the coefficients such 
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solutions can be taken so that the Jacobian of / evaluated at them is 
nonzero. 

By Theorem [6] this completes the proof of Theorem [2] 

3.3. Proofs of Corollaries [3] and 

Proof of Corollary 0 To prove Corollary [3] we follow the steps of the 
proof presented in subsection 13.11 however we highlight some differ¬ 
ences. We consider the change of coordinates given in (J9J) applied to 
system (J2D taking a 00k = b 00k = 0 for all k and instead of (ITU]) we obtain 

n 

9 = 1 + e ^ r l+j z k (a ijk cos* 9 sin 741 9 + bij k cos* +1 9 sin 7 9 ) , 

i-\-j-\-k =0 
n 

r = £ r'* +j, V' : (a ljk cos* +1 9 sin 7 9 + bij k cos* 9 sin 7 " 1 " 1 0) , 

i-\-j-\-k =0 
n 

Zi = £ ^2 c tijk r l+j z k cos* 9 sin- 7 9, 

i+j+k =0 

for i — 1, 2,..., d. 

As r does not appear in the denominator of 9, if £ is sufficiently 
small 9 > 0 for every ( 9 , r, z) in a ball B of an arbitrary given radius 
around the origin of W i+ ' 2 . Now in the ball B the variable r can be 
approximated to the zero as we want, this cannot occur working with 
the cylindrical annulus A of subsection 13.11 From here the calculations 
are done analogously to section 13.11 and we obtain the same maximum 
number of zeros of the averaging function for the new system (ITT]) with 
o-ook = book = 0. □ 

Proof of Corollary [J} The same above argument can be used for prov¬ 
ing Corollary [I] and we follow the steps of the proof presented in sub¬ 
section More precisely, we apply the change of coordinates (JUJ) to 
system (JUJ) considering a 0ok = a 00k = b 0 o k = /?oo k — 0 for all k in order 
to obtain an expression for 9 in both systems y > 0 and y < 0 with 
denominator that does not depend on r. So for £ sufficiently small we 
have 9 > 0 for every (r, 9, z) in the ball B of the proof of Corollary [HJ 

Since the exponent i + j of the variable r on the polynomial /i, given 
in (fl7|) . is at least one because i is odd, we have 


(18) 


/i(r, z) = Air + A 2 r 2 + ... + A n r n 
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n—p 

'y ' y ' z ( a ijk + ( —1 yoiijk) I(i+l,j) 

k=0 i + j = P 
i odd 

n—p 

+ Z k {pijk — ( — ^ypijk) 

k=o i-\- j = p 
i even 

for p = 1,2,... ,n. 

Consequently we obtain f\ = rfi where 

A O', z) = A 1 +A 2 r + ...+ A n r n_1 . 

As r > 0, to know the solutions of (/i, / 2 ,..., fd+i) = (0,..., 0) is 
equivalent to solve (A, h ,..., fd+i) = (0, But, A is a poly¬ 

nomial on variables r and z of degree n — 1 and the functions A +i 
given in (H71) are polynomials on variables r and z of degree n for each 
t — 1, 2,..., d. By Bezout’s theorem / has at most n d {n — 1) solutions 
(with r > 0). 

To prove that this number is reached we choose a particular example. 
This step of this proof is done as in subsection 13.21 for obtaining the 
maximum number of zeros of the averaging function for the new system 
(JT5]) with a 0 ofc = «oo k = b 00k = A) 0 fc = 0. Again in the ball B the 
variable r can be approximated to zero as we want. 

□ 


with 

A p = 
(19) 
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